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ABSTRACT 

The  effects  of  boundary  perturbations  on  finite-amplitude 
acoustical  standing  waves  in  rectangular,  rigid-walled 
cavities  were  investigated  using  non-linear  theory.   When  a 
high  amplitude  standing  wave  of  frequency  ca  is  generated  in 
a  cavity,  non-linear  effects  v;ill  cause  a  stimulation  of 
certain  normal  modes  whose  resonance  frequencies  are  integer 
multiples  of  oo.   Previous  experimental  observations  revealed 
that  there  could  be  excitation  of  other  normal  modes,  not 
belonging  to  the  family  of  the  driven  mode,  which  v/as  not 
predicted  by  the  non-linear  theory. 

The  purpose  of  this  research  v/as  to  investigate  the 
possibility  that  deviations  from  the  idealized  geometry  could 
account  for  these  observations.   Of  the  various  mechanisms 
possible,  this  work  investigated  the  possibility  of  these 
unpredicted  excitations  occurring  through  a  non-linear 
mechanism. 

The  standing  waves  that  exist  in  an  ideal  cavity  must  be 
corrected  when  the  boundaries  are  irregular.   The  non-linear 
interaction  between  these  standing  waves  and  the  corrections 
was  studied.   The  ability  of  this  interaction  to  excite 
standing  v/aves  other  than  those  predicted  in  the  ideal  case 
was  verified.   A  specific  example  was  worked  out  demonstrating 
an  unpredicted  excitation,  the  strength  of  v/hich  was  on  the 
order  of  the  magnitude  of  the  boundary  perturbation  parameter. 
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I.   INTRODUCTION 

The  purpose  of  this  research  was  to  investigate  some 
of  the  effects  of  boundary  wall  perturbations  on  finite- 
amplitude  standing  waves  in  a  rectangular  cavity.   The 
investigation  was  prompted  by  an  examination  of  the  experi- 
mental results  of  Coppens  and  Sanders  [1]  and  the  research 
of  DeVall  [2] ,  which  suggested  the  existence  of  the 
excitation  of  modes  other  than  those  belonging  to  the  family 
of  the  driven  mode. 


II.   BACKGROUND 

An  acoustical  wave  traveling  in  a  medium  will  have  its 
form  distorted  as  predicted  by  the  non-linear  hydrodynamic 
equations.   If  the  medium  is  absorptive  then  only  waves  of 
relatively  high  amplitude  will  develop  distortion  before 
becoming  completely  attenuated  due  to  the  absorptive  process. 

At  the  United  States  Naval  Postgraduate  School,  Coppens 
and  Sanders  [1],  DeVall  [2],  and  others  [3,4]  have  dealt 
with  the  study  of  finite  amplitude  waves  in  rigid-walled 
rectangular  cavities  and  tubes.   One  interesting  result  of 
these  cavity  experiments  was  the  appearance  of  excitations 
of  modes  which  were  not  family  members  of  the  driven  mode. 
For  example,  assume  that  the  cavity  is  driven  in  the  (0,1,0; 1) 
mode,  which  indicates  the  excitation  is  of  the  form 
cos  k  y  cos  cot  ,   where 


y   L  _ 


and 


Y  c 


Then  a  non-familial  mode  is  any  which  does  not  have  the 

form   (0,n,0;n)  .   It  was  proposed  [5]  that  the  presence 

of  these  non-familial  modes  could  be  attributed  to  the  cavity 
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walls  being  non-ideal  (not  perfectly  flat)  and  that  perhaps 
the  previous  small  disparity  between  theoretical  predictions 
and  experimental  observations  might  be  partly  explained  by 
the  presence  of  these  non-familial  modes. 

The  non-linear  wave  equation  applicable  to  this  problem 
can  be  expressed  [6]  as 


r^  2_^2     P    _     1  9^   r   /   P   ^2  ^  /  0,2, 
^0  Pl  2   -     2  ~2  ^"Y  (7 2^   "^  ^SZ^     ^ 


Po^o"     "  9^"    '0^0^     ""o 


.  1  ^  2  „2  r  /   P  ,2         ,^0,  2, 

+  2  ^0   ^   f( 2^   "  ^^^  ^   • 


(2.1) 


2 


is  defined  as 


ul  =  d'  -  ^ 


where  L      represents  the  loss  terms  appropriate  to  each 
frequency  component  in   P  .   The  RHS  of  the  non-linear 
equation  may  be  interpreted  as  a  forcing  function  (consisting 
of  both  classical  and  non-linearly  generated  terms)  which, 
in  first-order  perturbation  theory,  is  the  classical  solution 
for  the  pressure.   The  first-order  perturbation  solution  may 
be  considered  to  describe  the  effects  of  self-interaction  of 
the  classical  solution,  whereas  higher-order  perturbation 
solutions  describe  the  effects  of  interaction  of  the 
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of  the  non-linear  solution  with  itself.   If  the  cavity  is 
driven  at  frequency   co  ,   non-linear  terms  force  the 
existence  of  all  terms  containing  the  form  nco  .   Each  such 
term  whose  frequency  is  close  to  the  resonance  frequency  of 
a  standing  wave  of  the  cavity  and  whose  spatial  behavior 
matches  that  standing  wave  can  strongly  excite  that  standing 
wave  to  a  degree  determined  by   Q  ,   the  quality  factor  of 
the  resonance,  and  by  the  difference  between   nw   and  the 
resonance  frequency  of  the  standing  wave.   In  a  rectangular 
cavity,  it  can  be  seen  [1]  that  if  the  classical  solution 
is  the   (m  ,m  ,m  ;q)   mode,  then  the  non-linear  wave  equation 
will  force  the  strong  excitation  of  all  modes 
(nm  ,nm  ,nm  ;nq)   which  constitute  the  family.   Any  mode  of 

X     j/^     z  ■  ■ 

frequency  noo  which  is  not  a  member  of  the  family  should 
not  be  strongly  excited.   The  purpose  of  this  research  is  to 
see  if  the  presence  of  wall  irregularites  can  explain  how 
non-family  members  may  be  strongly  excited,  and  to  present 
an  example  to  support  this  theory. 
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III.   THE  NON-LINEAR  PROBLEM 

A.   CAVITY  DESCRIPTION 

Assume  a  perfectly  rigid-walled  rectangular  cavity  has 
one  wall  perturbed  such  that  the  cavity  dimensions  are 
L  [l+6f(y)],  L  ,  L   .   A  top  view  depicting  the  XY  plane 
and  showing  the  boundary  irregularity  is  given  in  Fig.  1 
below.   Also  assume  the  perturbation  on  the  boundary  is 
small  compared  to  the  cavity  dimensions,   l6f (y)  j  <<  1  . 


FIGURE  1 


The  cavity  is  to  be  excited  by  a  source  near  one  of  the 
corners  in  such  a  way  that  the   (n,m,£;q)   mode  is  excited. 
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B.   CLASSICAL  SOLUTION  FOR  IDEAL  BOUNDARIES 

For  a  rigid-walled  rectangular  cavity  with  ideal 
boundaries  (6=0),  the  classical  linear  solution  for  the 
pressure  p^   is  subject  to  the  following  conditions, 
namely 


a\  =  0  (3.1) 


and      ^Pq'^O   "   °     ^^   ^"^'^x  (3.2) 

y=0,L 
z=0,L^   , 


where   n_   is  the  local  normal  to  the  ideal  surface.   The 
solution  for   p-^   in  terms  of  Mach  number  is 


Pn 

TT     =   M  cos  k  X  cos  k  y  cos  k  z  cos((jo  t+c}))   ,   (3.3) 

z  X       y       z       cj 


where 


k    =   ?^  (3.4) 

X         Li 

X 


miT 

S    =  IT 


y 


k   =  ■= — 
z      L 
z 
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and 


(^)2   =   k2   =.   /k  2  +  k  2  +  k  2     ^  (3  5j 


where   w    is  the  frequency  associated  with  the   (n,m,£;q) 
mode. 


C.   CLASSICAL  SOLUTION  FOR  REAL  BOUNDARIES 

Since  the  cavity  is  assumed  to  have  one  perturbed  wall, 
it  is  reasonable  to  expect  the  solution  for  real  boundaries 
to  differ  from  that  for  ideal  boundaries.   The  classical 
linear  solution  for  real  boundaries,   P^  ,   may  then  be 
considered  a  summation  of  the  classical  linear  solution  for 
ideal  boundaries  plus  perturbation  correction  terms  due  to 
the  irregular  boundary: 


Pq   =  Pq  +  -SPi  +  <5  P2  +  .  .  . 


As  long  as  the  magnitude  of  the  boundary  perturbation  is  kept 
small,  terms  higher  than  first-order  can  be  considered 
insignificant,  so  that 

P.   =  Po  +  5p,     (to  first  order)         (3.6) 

and  must  be  a  solution  to 


□\  =  °  '^•'' 


15 


with  boundary  conditions 

VPQ-n   =   0     at   x=0,  L^[l  +  6f(y)]  (3.8) 

y=0,  L 

y 

z=0,  L     , 

where   n  ,   the  local  normal  to  the  real  surface,  is  obtained 
by  taking  the  divergence  of  the  equation  for  the  perturbed 
boundary,  i.e.,   V • {x  -  L  [l+6f (y) ] }  .   Thus,   to  first  order 
in   6, 


^   •  ^       T  r  3f (y)  ^ 


and  when  the  above  is  used  in  (3.8)  the  result  is 


8P 


0 


8x 


9P, 


-  L  6  — ^  9f  (y) 

x=L  [l+6f(y)]     "^   ^^    ^^ 


=   0 
x=L  [l+6f (y)] 

(3.9) 


D.   THE  APPLICABLE  NON-LINEAR  EQUATION 

A  real  wave  in  a  real  cavity  will  of  course  have  non- 
linearities  introduced  as  it  propagates  through  the  cavity. 
Let   P,   be  the  correction  term  to  the  pressure  due  to  the 
non-linearities  of  the  medium.   The  total  pressure   P   is  tlien 


P   =   pQ  "^  ^1  (3.10) 
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When  (3.10)  is  inserted  in  (2.1),  the  LHS  of  the  non-linear 
equation  becomes 


P  P 

2_2     "^0   .    2^2    ^1 


^0  Ql   2  "■  ^0  Dl 


The  first  term  above  is  known  to  be  zero  from  (3.7),  so  the 
LHS  of  (2.1)  becomes  simply 


2^2    ^1 
^0  Dl 2 


Now  if  standing  waves  are  to  be  produced  in  the  cavity,   P^ 

must  consist  of  terms  of  the  form  cos  k  x  cos  k  y  cos  k  z  cos  ot 

X       y-*       z 

Using  the  basic  definitions 

^0     "Pq  dt 


and 


3o  =  '*o 


requires  that  U^   be  of  the  form   sin  k  x  sin  k  y  sin  k  z  sin  wi 

(J  X       y       z 

When  this  information  is  used,  the  RHS  of  (2.1)  will  consist 
of  four  terms : 


3^  2 

1.   -  •2— :r-  (cos  k  X  COS  k  y  cos  k  z  cos  cot) 


17 


92  2 

2.   -  — y  (sin  k  X  sin  k  y  sin  k  z  sin  wt) 
3t  ^       ^ 


2  2 

3.   V   (cos  k  X  cos  k  y  cos  k  z  cos  wt) 
X       y-'       z 


2  2 

4.   V   (sin  k  X  sin  k  y  sin  k  z  sin  wt) 


Although  the  expansion  of  these  will  produce  many  terras,  it 
must  be  noted  that  only  some  will  produce  resonant  standing 
waves.   All  the  other  terms  may  therefore  be  considered  to 
give  a  small  contribution  when  compared  to  the  resonant  terms , 
and  thus  may  be  neglected.   This  fact,  when  used  with  the 
trigonometric  identities 


2     1 
cos  a  =  2"  [1  +  cos  2  a] 


and 


2     1 

sin  a  =  J  [1  -  cos  2a] 


reduces  the  RHS  of  (2.1)  to  a  single  term,  and  after  a  good 
deal  of  manipulation  the  non-linear  equation  may  be 
approximated  by 


P  2      P 

^1    ^   _  y  +  1  8l_  (__J 

2   "      2     2 


^22    1   ^  _  y  + 1  a   / o_.  2       ,3  ..X 


The  governing  non-linear  equation  is  now  in  a  relatively 
simple  form,  and  Appendix  I  develops  the  general  form  of 
the  solution  for   P-,  . 
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IV.   THE  PERTURBED  BOUIJDARY 

A.   DEVELOPMENT  OF  A  BOUNDARY  CONDITION  ON  p. 

Recall  that   p,   is  the  first-order  correction  term  due 

to  one  boundary  being  irregular.   Although  it  is  known  that 

for  a  perfectly  rigid-walled  cavity  the  pressure  must  be  a 

maximum  at  the  walls,  this  fact  may  not  be  conveniently 

used  with  one  wall  having  an  irregular  form  such  as 

L  [l+6f(y)]  .   It  would  be  highly  desirable,  therefore, 

to  develop  a  boundary  condition  on   p,   to  be  applied  at   L   . 

This  may  be  accomplished  as  follows: 

Using  the  incremental  form  for  a  Taylor  series  expansion 

on   Pp.   evaluated  at  the  real  boundary   L  [l  +  6f(y)]   produces 
u  X 


x=L  [l+5f (y)] 


3P, 


+ 


8x 


x=L 


a^p 

1  '^   0 


x=L 


[L^,6f(y)]"  + 


x=L. 


X 


Substitution  of 


Pq  =  Po  -^  "^Pl 


(3.6) 


on  the  RHS,  and  taking  the  partial  derivative  with  respect  to 
X  on  both  sides  gives 
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Hi 

dx 


x=L^[l+6f (y)  ] 


9P, 
dx 


+   6 


x=L 


8p- 
9x" 


X 


x=L. 


8  p. 


9x 


x=L 


Lj^6f(y)  +  ... 


.  (4.1) 


X 


Following  exactly  the  same  procedure  for  the  y   coordinate 
gives 


ap< 


dy 


3P, 


y=L^[l+6f (y)] 


3y 


+  6 


x=L 


8p. 
W 


x=L 


X 


9  P, 


+ 


By' 


x=L. 


L^6f(y)  +  ... 


X 


.  (4.2) 


Recall  that  (3.9)  required 


8P, 


dx 


x=L  [l+6f (y)] 


T  R   ^^0  9f (y) 


=  0 


x=L  [l+6f (y)] 


When  (4.1)  and  (4.2)  are  used  in  (3.9)  the  result  is 
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9P. 


8x 


68P- 


x=L. 


dx 


9  P< 


x=L^    9x 


x=L 


L^6f(y)  + 


L  6|iiXL[!^ 
X  8y    ^8y 


+ 


X=L 


6  8p. 

9y~ 


9  Pq 


x=L     8y 


X 


x=L 


L  6f  (y)  + 


]=0 


Collecting  like-order  terms  in  6  then  gives 


9P- 


9x 


+ 


a^Po 


x=L     8x' 

X 


X=L 


L  f(y)  -  Liilil!^ 

X  ^^'  X  8y    3y 


=  0, 


x=L. 


or,  as  a  boundary  condition  on  p,  applied  at  L  , 


dp. 


dx 


9  Po 

=  L  [-f (y)  — ^ 
x=L      ^        8x^ 


3f  (y)  ^Ppj    , 
x=L    ^i-    ^^     x=L 

X  X 


(4.3) 


B.   A  FOURIER  SERIES  REPRESENTATION 

It  will  be  desirable  for  further  development  to  represent 
the  RHS  of  (4.3)  as  a  Fourier  series  in  cosineS/  so  that  p, 
may  be  expressed  as  a  summation  of  normal  modes.   First, 
rewrite  (4.3)  as 


A1L^[-f  (y) 


2 

9  PO 


8x 


x=L 


+  9f (y)  ^^0 


ay  ay 


x=L 


].   (4.4) 


x=L 


X 
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Let  the  RHS  of  (4.4)  equal  Scos  k  z  cos  cot  and  express  it 
as  an  infinite  Fourier  series  as  follows: 


A  00 

S(y)cos  k^z  cos  cot  =  [-^  +  2  A^  cos  ^^]  cos  k^z  cos  cot, 

N=l        y 


where  the  A  are  given  by 

L 

Aj^  =  j7-   /   S(y)  cos  ^  dy,      N  =  0,1,2, 

y  0  y 


Then  (4.4)  becomes 


N=l  y 

x=L 

X 


3x^~^^ 
P  c 
0  0 


The  solution  for  p^  is 


p,         cos  k  'x 

— J  =  -  E  A^  cos  -~-  cos  k  z  cos  ojt    (4.5) 


DC        N  k  'sin  k  'L  y 

^0  0  X      XX  ^ 


where  for  N  =  0  the  coefficient  is  A„/2,  and  k  '  is  such  that 

0   '      X 


X      7  c       L^       z 


Note  that  k  '  is  not  in  the  form  of  a  standing  wave  in  the 
X  ^ 


cavity  required  by 
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X 


Thus  the  terms  involving  k^'  in  (4.5)  must  also  be  expressed 


X 

as  a  Fourier  series: 


cos  k  'x        B      00 

X-  +   Z   B   cos  — —  ,  (4.6) 

k^'  sin  \'    \  ^         M=l   ^      ^x 


where  the  B   are  given  by 


M 


^     X   cos  k  'x        -. 

B»   =   j^   / _  cos  ^  dx,   M  =  0,1,2, 

X  0     k  '  sin  k  'L        x 

X  XX 

When  the  results  of  (4.6)  are  inserted  into  (4.5),  p,  is 
expressed  as  a  double  infinite  series: 


p 


^  =  -  Z  E  B   cos  — y—  A^  cos  —Y~   cos  k  z  cos  cot,  (4.7) 


PqC^       N  M  X  y 


where  again  each  of  the  zero  coefficients  (A  ,  B  )  are  to  be 
^  0    0 

divided  by  a  factor  of  two. 

The  formalism  has  now  been  developed  whereby  given  any 
driving  mode  and  boundary  irregularity,  p,  -  the  first-order 
perturbation  correction  term  -  may  be  expressed  in  terms  of 
the  spatial  components  of  the  normal  modes  of  the  ideal 
cavity. 
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V.   A  SPECIFIC  EXAMPLE 

A.   THE  COS INUSOI DALLY  PERTURBED  BOUNDARY 

As  a  specific  example  utilizing  the  theory  developed 
in  the  preceeding  sections,  assume  the  cavity  is  driven  in 
the   (0,1,0; 1)   mode.   Thus 


Pn^ 


M  cos  k  y  cos  cot 


(5.1) 


0^0 


Further  assume  that  the  boundary  is  perturbed  cosinusoidally 
such  that 


f  (y) 


7T 

cos  ~—  y 


(5.2) 


The  boundary  condition  on  p-,   then  becomes 


9_  (_ix.) 

3^   p  c  ^ 
^0^0 


x=L 


M  T-^  TT  k   sin  ~-  sin  k  y  cos  wt 
L     y     L       y-' 
y    ^      Y 


T^ 


which,  when  used  with  the  requirement  that  k  ,  =  r—  , 

y    J-i 

^        y 

becomes 


^    (— X) 

^0^0 


x=L 


X 


M  _2i  I  k   cos  wt  -  M  =^  i  k   cos  ^^  cos  ^t 
Ly  2  y  Ly  2  y     Ly 
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By  inspection,  the  solution  for   p,   must  be 


—^-      "   ^   —  2    k^  sin   k    'L    ^°^    ^t 


pQ^o  y        X  XX 


L  k       cos    k    "x  ^ 

o-i/XTry  X  2  Try  ^  /,--^x 

y         X  XX  y 


where   k  '   satisfies  the  equation 


K-=1  (5.4) 


and  k  "  satisfies 


y 


Solution  for   p-.   allows  the  creation  of   P^  =  p^  +  6p,   , 


P«  L     k   cos  k  'x 

r-  =   M  cos  k  y  cos  cot  -  6M  ~ —  -^  r—\ = i — rr-   cos  cot 

^  _  2  y^  L   2  k  '  sin  k  'L 

PqCq  -^  y    X       X  X 


L     k  cos  k„"x  o 

,r.,/XTTy       X  2  Try  .           it-    r\ 

^  ^^^  L"  2  k^  sin  k  "  ^°^  L   ^°^  ^^   •        ^^-^^ 

y    X      X  y 


Let  the  cavity  dimensions  be  designed  such  that  the   (0,2,0; 2) 
and  the   (1,1,0; 2)   modes  will  be  degenerate.   This  requires 
that 
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y 


and 


/_n_\  2    r_IL^  2   _   ,2a)^2 
X       y         ^ 


Solution  of  the  above  equations  yields 


j^  =   /3  (5.7) 

X 


for  the  desired  degenerate  modes.   When  (5.7)  is  inserted 
into  (5.6),  the  result  is 


^0  1      TT     ^v     """^    ^x'^ 

75-     =      M   cos    k  y    cos    cot   -    6M  -^  r— ^ < r: — rr-  cos    wt 


PqCq  ^  /3"  X  X     X 


-  k  cos   k    "x  o 

.    X  II      1      f^       Y                 ^  2Try  .  /con 

•^   ^^  7=-  2  k^  sin   k    "L    ^^^  L       ^^^    ^'^       •               ^^-^^ 

/3  X                XX  y 


This   must  be   squared  when  substituted  into    (3.11). 
Recalling   that 

Pq    =   Po    +    «Pl       '  '3.6) 


26 


then  the  non-linear  term  has  the  form 


^0      ^               Pfl   ^        Pfl      P]       2    P] 
( ^)   =  ( 0^)   +  2«( Oj)  ( i^)  +  «2( 1^, 

"o^a     pQ^o       Po'^o   "o'^o      "o'^o  ■ 


2 

The  third  term  on  the  RHS  above  is  of  order  6   and  may 

therefore  be  neglected.   The  first  term  becomes 


( ^)  =  M   cos   k  y  cos   cot. 

2  y-' 

DC  -^ 

^0  0 


Using  the  trigonometric  identity 


2      1 
cos   a  =  2-  [1  +  cos  2a] 


2 
produces  the  following  for  p   : 

(- ^)       =  ilL_  [1  +  cos  2k  y]  [1  +  cos  2cot]  . 

^0   0 

Since  the  only  significant  contribution  will  come  from 
terms  near  resonance,  those  terms  not  involving  a  product 
of  cos  2k  y  cos  2a)t  may  be  neglected.   Thus 

Pn   2    .,2 
( ~)       =  ~   cos  2k  y  cos  2a)t   .  (5.10) 

^0  0 
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similarly,    keeping  just   those   terms  which  will    lead   to 
resonance  produces 

Pa           Pi  J.1         o          o  1           ^        cos   k    'x 

ox           0              1      _  Y+1    XI/2/0    \^  1      IT      y                  X                    ,               ^    ^ 

26   2"  J  ~    ~   "2~  ^       ^^'    —  ?" ^^^    ^   ycos2cot 


»,jli  o  o        t        _'K  COS/C        X 

+  -^^5—  6m       (2co)       ■ —  -r-  rf-ir  cos    k   y    cos    2a)t     . 

^  /I  ^    ^x      sin   k    "L  ^ 

^      ^  (5.11) 


Inserting   the   results   of    (5.9    -    5.11)    into    (3.11)    gives 


c 


2_^2  1  Y+1      2      2  ^,  _    ^ 

0.    FIl  7  "^      2~   '       ^      ^°^    ^^v^   ^ 

PO^O 


,     ^  r^  «  -|  JC  COS         JC  X 

-   ^  6    M      (203)      ~  ~  -^ ~ COS   k  y   cos    2'a.it 

/3   ^   k    '    sin  k    'L  ^ 

X  XX 

,  p  _     ,  k        cos   k    "x 

+  ^V^  6    M      (2w)      —  ^  -^ cos   k  y   cos    2ujt. 

/J       k    "    sin  k    "L  ^ 

^  (5.12) 


The   terms    cos   k    'x   and   cos   k    "x  must  be   expressed   as    a   cosine 

XX  ^ 

series  of  the  form  given  by  (4.7) .   When  this  is  done,  and 
only  those  components  of  argument  ttx/L   are  retained,  then 
all  terms  driving  the  (1,1,0; 2)  mode  will  be  isolated,  as 
desired.   The  results  are 


k^'   =   -rr/Ly  (5.13a) 
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COS  k  'x  =  0.535  cos  - —  ,  (5.13b) 

X  Li 

X 


sin  k^'L^  =  0.971  ,  (5.13c) 


^J'    =    i*^^  (5.13d) 

X      —J- 

y 

cos  k^"x  =  -3.677  cos  ^   ,  (5.13e) 

X  JLi 

X 


and 


sin  k  "L   =  123.1  (5.13f) 

X   X 


where  use  was  made  of  the  identities 


sin  iA  =  i  sinh  A 


and 


sinh  A  =  ^  [e   -  e   ]  . 


When  the  results  of  (5.13)  are  inserted  into  (5.12),  the 
result  is 


P 

C-    n  T      T  -      o~  W      to      COS    2k   y   cos    2 cot 

0    LJi-"  ^2  y 

Po   0 
-    61.83   ^^—  M      CO      cos   J—  cos   ~  cos    2a)t,  (5.14) 

£.  Li  Li 

X  y 
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where  the  first  term  on  the  RHS  is  called  the  finite 

amplitude  term,  and  the  second  term  on  the  RHS  shows  the 

effects  of  the  boundary  irregularity  on  the  first  order 

non-linear  perturbation  correction. 


B.   THE  SOLUTION  FOR  P 

Appendix  I  shows  the  general  solution  for  a  function 

2    2 

operated  on  by  c   Q^,  •   Using  those  results  /  the  solution 

for  P,  for  this  example  becomes  immediate:   for  the  finite 
amplitude  term,  designated  by  subscript  t, 


^1  Y+1  ^f  ^^"  ^f  2   2 

( K)       =  ^^i  — 7f-  ir    00   cos  2k  y  cos  (2ajt+3.)  ,  (5.15) 

^2  2      /o  \ 2  y            f 

PqCq   f  (2a))  ^ 

and  for  the  boundary  irregularity  induced  term,  designated 
by  subscript  p. 


Pt  y-t   9   9  Q^  sin  3^ 

( i^)   =  -1.83  6  -i^  M^  03  -E ^  cosk  xcosk  ycos  (2a3t+3„) 

PqCq   P  (20))''       X     y  P 

(5.16) 


When  (5.15)  and  (5.16)  are  combined,  they  form  P, ,  the  first 

order  non-linear  perturbation  correction  to  the  pressure. 

P,   is  a  non-linearly  generated  term  which  results  from 

■^f 
the  existance  of  a  010  standing  v/ave  in  the  cavity.   P, 

P 
exists  because  the  irregular  boundary  forces  a  non-linear 

interaction  between  the  correction  term  p,  and  the  010  mode. 

It  should  be  noted  that  there  exists  an  additional 

correction  term  to  the  total  pressure.   This  is  the  classical 


30 


correction  to  the  020  standing  wave  which  exists  in  a 
cavity  with  an  irregular  wall.   This  additional  correction 

term  will  be  on  the  same  order  of  magnitude  as  P,   [7]  , 

P 
and  must  be  taken  into  account  when  finding  the  total  effect 

of  the  irregular  boundary  on  the  pressure.   The  following 

discussion  of  results  will  consider  only  the  correction  term 

generated  by  this  example. 

C.   RESULTS 

In  order  to  best  illustrate  the  magnitude  of  the  effect 
of  the  non-familial  mode  in  this  example,  recall  that  from 
(5.15)  and  (5.16) , 


P  ,    Q      sin    3^r      o      o 

1  Y+i  f  '^f       ,,2  2  ^,  /O      a.1   «      \ 

2"  ^      2 2 — ■   ^      ^      ^°^         y^    ^°^     (2Lot+3^) 


Po^o  ^^'^^ 


^^   Q      sin    3         2      2 

-    1.83  5  -^ ^ 9—^  M      0)      cos   k  x  cos    k  y   cos    (2cot+3 

^     (2a))  ^  ^       ^  ^ 


Dividing  both  sides  of  the  equation  by 


y+1  M^ 
2   4 


allowing  the  quality  factors  Q^  and  Q   to  be  equal,  and 
assuming  measurements  are  made  at  one  of  the  corners  so  that 


cos  k  X  =   cos  k  y   =   cos  2k  y   =   1, 
X  y-^  y 
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gives 


_i Q_Q 


^-Li =   sin    6^   cos    (2oot+3£)    -    1.836sin   B      cos    (2cot+B    ) 


Y+1  M 


Let  the  LHS  be  designated  by  C.   Additionally,  introduce 


H   =   2    ^   P   (-i E) 

Q^  +  Q     w 


as  a  measure  of  the  overlap  of  the  curves  obtained  by 
plotting  the  magnitude  of  C  versus  a  function  of  frequency 
for  two  cases: 


1  -   6=0,   corresponding  to  no  non-familial  mode 

existance,  and 

2  -   6  7^  0 ,   corresponding  to  a  degeneracy  between  the 

(0,2,0; 2)  and  (1,1,0; 2)  modes. 


These  curves.  Figs.  2-6,  allow  a  comparison  of  the  magnitude 
of  the  P,  term  with  and  without  the  presence  of  a  non- 
familial  mode.   The  function  of  frequency  on  the  abscissa 
is 


F  =  cot  Bf/ 


where  B  is  as  defined  by  (I.l).   Plots  are  provided  for  five 
different  values  of  H.   Additionally,  for  the  case  H  =  -1 
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there  is  a  plot  labelled  H_  =  -1  which  shows  the  result  of 
selecting  a  different  corner  for  measurements,  one  180° 
out  of  phase  with  the  plot  labelled  H   =  -1.   6  was  arbi- 
trarily chosen  as  0.1. 

From  the  graphs  it  is  clear  that  the  presence  of  a  non- 
familial  mode  does  affect  the  total  solution  for  the  pressure 
by  an  amount  proportional  to  6,    the  magnitude  of  the  boundary 
perturbation . 
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VI.   CONCLUSIONS 

Non-linear  theory  has  been  applied  to  standing  waves 
in  a  rigid-walled  rectangular  cavity  with  a  perturbed 
boundary  in  an  attempt  to  find  one  possible  mechanism  for 
the  excitation  of  resonant  modes  other  than  those  belonging 
to  the  family  of  the  driven  mode.   It  was  observed  that  such 
an  excitation  does  exist  if  the  boundary  perturbation  and 
the  dimensions  of  the  cavity  are  favorably  chosen.   If  the 
fractional  perturbation  of  the  boundary  is  given  by  6 ,  then 
the  correction  terms  are  of  that  same  order. 
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APPENDIX  I 

It  IS  desirable  to  show  that  when   c   Q^    operates 
on  a  function  of  the  form  cos  k^x  cos  k^^y  cos  k^z  cos  (wt+c}))  , 
the  result  will  be  of  the  form 


%■  sin  3  ^°^  ^x^  *^°^  \^   ^°^  ^z^  cos(ajt+(|)-3)  , 
where   3   is  an  angle  defined  such  that 


c  k 
I  /    [1   -    (-^)^]   =   tan  3   .  (I.l) 


Begin  by   recalling   that,    by  definition. 


r.2 

22  2„2d  (jo8  -.t.«. 

°P  Dl      =    '^p  '    -  7^2  -  q  at     •  '^-^l 


Performing  this  operation  on  a  function  of  the  form 

cos  k  X  cos  k  y  cos  k  z  cos  (oot+d))   produces  three  terms,  all 
X       y-'       z  ^    ^ 

with  the  common  factors   cos  k  x  cos  k  y  cos  k  z  .   Canceling 

X       y-^       z  ^ 

these  like  terms  from  both  sides  of  (1.2),  and  using  the 
fact  that 


k^   =   k^+k^+k^ 
X     y     z 


requires  that 
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-c_    k      cos(aJt+(f))    +   oj^cos  (oot+(|>)    +   ^  sin(wt+4>) 


W^     1  ,   ^ 


or,  collection  terms, 


2  2 

2      "^n^  2  2    , 

CO  [1  -   P  ^  ]  cos(aJt+(|))  +  ^  sin(cot+({))  =  — J^-r-   cos  (cot+(t)-3) 

^3^^  Q         ^     Q  sin  P 

(1.3) 


Rearrangement  of  (I.l)  gives 


2  2 

c  ^y. 

[1  -   P  ^  ]   -    cos  B 


2  ■■      Q  sin  3   ' 

0) 


and  substitution  of  this  result  into  (1.3)  requires 


2  2  2 

^  fx|-|  cos(a)t+(^)  +  ^  sin(a3t+cf))  =  ^  ^j^  cos  (a>t+(j)-3)  . 


2 

Canceling  the  common  factors   oo  /Q  ,   and  multiplying  both 

sides  by   sin  3   then  requires  that 


cos  3cos(a3t+(}))  +  sin  3  sin((jot+({))  =  cos  (a)t+(})-3)   .      (1.4) 
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Using  the  trigonometric  identity 


cosA  cos  B  +  sin  A  sin  B   =   cos (A  -  B) 


immediately  shows  (1.4)  to  be  satisfied.   Thus  a  general 

2   2 
solution  to   c   Q    when  operating  on  a  function  of  a 


specific  form  has  been  shown: 


2   2 

c   pIt   cos  k  X  cos  k  y  cos  k  z  cos(cot+(})) 
p  UlL       X       V       z 


2    -, 

^^^r : ^  cos  k  X  COS  k  y  cos  k  z  cos  (ajt+(J)-3) 

Q  sm  3      X       y-^       z 


,       c  k  „ 
where   tan  3  =  ^  [1  -  (-^)  ]   • 
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